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Abstract 



In this paper we obtain a Feynman-Kac formula for the solution of a 
fractional stochastic heat equation driven by fractional noise. One of the 
main difficulties is to show the exponential integrability of some singular 
Mh . nonlinear functionals of symmetric stable Levy motion. This difficulty will 

(-H ' be overcome by a technique developed in the framework of large deviation. 

This Feynman-Kac formula is applied to obtain the Holder continuity and 
moment formula of the solution. 

^ 1 Introduction 

> ■ 

f>^ , Let < a < 2 and let A — J2'i=i '§^ be the Laplacian. In this paper, we shall 

r~|~ I obtain a Feynman-Kac type formula for the following stochastic equation driven 

"^ ' by fractional noise 

o . 




dtdxi- ■ ■ dxd (1-1) 



where W(t, x) is a fractional Brownian sheet with Hurst parameters Hq in time 
k><( ■ and {Hi , . . . , Hd) in space, respectively. More specifically, for the solution u(i, x) 

}_( , to the above equation, we can write down the following Feynman-Kac formula 



u{t,x) — E 



X 



/(XHexp / / 5{X^_,-y)W{dr,dy) 



(1.2) 



where E^ denotes the expectation with respect to the symmetric a-stable Levy 
motion Xf , and 5 denotes the Dirac delta function. 

When -— — — in (11.11) is replaced by a (deterministic) continu- 

atoxi ■ ■ ■ oxd 
ous function c(i,x), this is the classical Feynman-Kac formula, which has been 
widely studied (see [?])• 
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In our above stochastic case, if — (— A)"/^ in (|l.ip is replaced by the 
classical Laplacian A (the case a = 2), the corresponding Feynman-Kac formula 
was studied first in [TT] in the case that all Hurst parameters Hq ,Hi,--- , Hj, 
are greater than or equal to 1/2 and then in |^ in the case Hq < 1/2. In 
both papers, the symmetric a-stable Levy motion is replaced by the standard 
Brownian motion. 

In |11| , the main difficulties to overcome are the following: The first one 
is to show the exponential integrability of some functionals of the Brownian 
motion so that (|1.2p is well-defined. The second one is to show that (|1.2p is a 
solution to (|1.1[) . In [3], even the existence of the stochastic integral in (II. 2p 
becomes a big challenge and was dealt with great care. We expect the similar 
difficulties will appear in our current situation. 

We shall follow the approach of [11] and [9] to approximate the noise in 
p.ip by smooth ones. However, to show the exponential integrability of the 
stochastic integral in (jl.2p we shall use a technique developed in j3]. This will 
largely simplify our computation even in the standard Brownian motion case. 
We shall also use Malliavin calculus to show that (jl.2p is indeed a weak solution 
to (IlT|) . 

It is straightforward to extend our results in the following two directions. 
First, — (— A)"/^ in (|l.ip can be replaced by a generator of more general Levy 
process or more general Markov process. Secondly, W{t,x) in (|l.ip can be 
replaced by a more general Gaussian field. We restrict ourselves to the specific 
case of (II. ip is to make our presentation simple and to present our idea and 
approach in a clear way. 

The paper is organized as follows. In Section 2, we very briefiy present 
some preliminary material on Malliavin calculus and stable Levy motion that 
we need to fix some notations. In Section 3, we give a definition to the stochastic 
integral in (jl.2p . This is necessary since the integrand involves the Dirac delta 
function. We also show the exponential integrability of the mentioned stochastic 
integral. In Section 4, we use Malliavin calculus to prove that (II. 2p is a solution 
to (|l.ip in a weak sense. Section 5 is an application of Feynman-Kac formula 
(|1.2p . We prove that the solution u{t, x) given by (jl.ip has a Holder continuous 
version. Wc also find the Holder exponent. The main part of the paper assumes 
that the stochastic integral involved in (|l.ip is in the sense of Stratonovich. 
In Section 6, we also discussed the Feynman-Kac formula where the stochastic 
integral in (jl.ip is in the sense of Ito-Skorohod. We also obtain a Feynman-Kac 
formula to represent the Ito- Wiener chaos coefficients of the solution. A formula 
for the moments of the solution is also given. Finally in Section 7, we present 
some lemmas used in this paper. 



2 Preliminaries 

2.1 Fractional Brownian motion and Malliavin calculus 

Fix a vector of Hurst parameters H = {Ho,Hi, . . . ,Hci), where Hi e (^)l)- 
Suppose that W = {W{t,x),t > 0,x E M''} is a zero mean Gaussian random 
field with the covariance function 



E{W{t, x)W{s, y)) = Rh„ (s, t) [] Rh, ( 



where for any H e (0, 1) we denote by Rh{s, t), the covariance function of the 
fractional Brownian motion with Hurst parameter H, that is, 

i?H(.,t) = i(itr+i.r-it-.n. 

In other words, T4^ is a fractional Brownian sheet with Hurst parameters Ha in 
time variable and Hi in space variables, i — 1, . . . ,d. 

Denote by £ the linear span of the indicator functions of rectangles of the 
form {s,t] X {x,y] in R+ x R"^, where {x,y] = {xi,yi] x • • • x {xd,yd]- Consider 
in £ the inner product defined by 

d 
{1(0, s] X {0,x] , l{0,t] X {a,y] )h = Rho («, t) ]_ _[ Rh, (xi ,yi). 

i=l 

In the above formula, if Xi < we assume by convention that I{o,xi] = ~I[xi,o)- 
We denote by H the closure of £ with respect to this inner product. The 
mapping W : I(o,t]x(o,x] ~^ W(t,x) extends to a linear isometry between H and 
the Gaussian space spanned by W. We will denote this isometry by 

W{(f))^ I I <j){t,x)W{dt,dx), 
Jo Jm'' 

ii (p E H. Notice that if (p and ip are functions in £, then 

X / (f>{s,x)ij{t,y)\s~t\^""~^T]\x,~y,\^"^-^dsdtdxdy, (2.1) 

JR2 xR2d f^ 



i=l 



where an = Y[i=o ^ii'^^i ^ '^)- Furthermore, H contains the class of mea- 
surable functions (j) on M.+ x M'' such that 

/ \(f,{s,x)(f>{t,y)\\s-t\^"«-^T]\x,^y,\^"^-^dsdtdxdy<oo. (2.2) 

jR2_xR2d f-^ 

We will denote by D the derivative operator in the sense of Malliavin 
calculus. That is, if i^ is a smooth and cylindrical random variable of the form 

F = fiWi<Pi),...,W{(f>n)), 



(f)i G H, f £ Cff (R") (/ and all its partial derivatives have polynomial growth), 
then DF is the "H-valued random variable defined by 

The operator D is closable from L'^{i}) into L'^{^; Ti.) and we define the Sobolev 
space D"'^'^ as the closure of the space of smooth and cylindrical random variables 
under the norm 



\\DF\\,^^ = ^E{F^)+E{\\DF\\l). 

We denote by 5 the adjoint of the derivative operator, determined by duality 
formula 

E{5{u)F)^E{{DF,u)^), (2.3) 

for any F e B)^'^ and any element u G L^(yi\ H) in the domain of J. The operator 
S is also called the Skorohod integral because in the case of the Brownian motion 
it coincides with an extension of the Ito integral introduced by Skorohod. We 
refer to Nualart [T^ for a detailed account on the Malliavin calculus with respect 
to a Gaussian process. If DF and u are almost surely measurable functions on 
R_(_ X R'^ verifying condition (j2.2p . then the duality formula (|2.3p can be written 
using the expression of the inner product in H given in ()2.1|) : 

E{S{u)F) ^an 

xE ( [ D,,,Fu{t,y)\s ~ t\''"'>-^f\ \x^ - y.l'^^'-^dsdtdxdy] . 

We recall the following formula, which will be used in the paper 

FWi(t>)=SiF(t)) + {DF,(t>)^, (2.4) 

for any (j) £% and any random variable F in the Sobolev space D^'^. 

2.2 Symmetric a-stable Levy motion 

In this section we recall symmetric stable distribution and symmetric a-stable 
Levy motion. For more general and detailed result about stable processes, we 
refer to P) . 

A random variable X is said to be symmetric a-stable if there is parameters 
< a < 2 and ct > such that its characteristic function 

and we will denote that X ^ S{a, cr). 

Notice that when a = 2, X is a Gaussian random variable. When a G (0, 2), we 

have E\X\P < oo if -1 <p < a and E\X\p = oo if p > a. 

A stochastic process {X{t),t > 0} is called symmetric a-stable Levy motion if 



(1) X{0) =0 a.s.. 

(2) X has independent increments. 

(3) X(t) — X{s) ^ S{a,{t — s)~) for any < s < t < oo and for some 
<a < 2. 

Throughout the paper C will denote a positive constant which may vary 
from one formula to another one. 

3 Definition and exponential integrability of the 
generalized stochastic convolution 

For any e > we denote by Pe{x) the d-dimensional heat kernel 
On the other hand, for any <5 > we define the function 

ips{x) == -I[o,5]{x). 

Then, (ps{t)pg{x) provides an approximation of the Dirac delta function d{t,x) 
as e and S tend to zero. We denote by W^'^ the approximation of the fractional 
Brownian sheet W{t,x) defined by 

W''\t,x)^ f f Mt-s)Pe{x-y)W{s,y)dsdy. (3.1) 

Jo Jr<' 

Fix a: G R'^ and t > 0. Suppose that X = {Xt, i > 0} is a d-dimensional 
symmetric a-stable Levy motion independent of W. We denote by Xf = Xt+x 
the symmetric a-stable Levy motion starting at the point x. We are going to 
define the random variable L L^ (5(Xf_^ — y)W{dr,dy) by approximating the 
Dirac delta function (5(Xf_^ — y) by 

Aliir, y)= [ ^s{t-3- r)p,{X^ - y)ds. (3.2) 

Jo 

We will show that for any e > and (5 > the function A^''^ belongs to the 
space "H almost surely, and the family of random variables 

V!j!= f j Ali{r,y)W{dr,dy). (3.3) 

converges in L^ as e and S tend to zero. 

The specific approximation chosen here will allow us in Section 4 to 
construct an approximate Feynman-Kac formula with the random potential 



W'^'^{t, x) given in (|4.1I) . Moreover, this approximation has the useful properties 
proved in Lemmas 17.41 and 17.51 We may use other types of approximation 
schemes with similar results. Also, we can restrict ourselves to the special case 
6 = e, but the slightly more general case considered here does not need any 
additional effort. 

Along the paper we denote by E^{'^{X,W)) (resp. by E^{<^{X,W))) 
the expectation of a functional $(X, W) with respect to X (resp. with respect 
to W). We will use E for the composition E-^ E^ , and also in case of a random 
variable depending only on X or VF. 

Theorem 3.1 Suppose that 2Ho + ^ J^t^iC^^, - 2) > 1. Then, for any e > 
and S > 0, A^'^ defined in iS. ^|) belongs to "H and the family of random variables 
V/^ defined in iS. 3\) converges in L^ to a limit denoted by 

Vt,,^ f f 5{Xf_,~y)W{dr,dy). (3.4) 

Conditional to X , Vt^x is a Gaussian random variable with mean and variance 

Var^(F,,,) ^c^H f f \r - s\^"»-^f\\x:. - Xl\^"''^ drds . (3.5) 

Jo Jo ,^1 

Proof Fix e, e', S and S' > 0. Let us compute the inner product 

(At:i,4f) ^ an I I p,{X^ - y)pAX^ - z) 
xips{t — s — u)ips' {t — r — v) 

d 

X |m - v\^"°~^ J]^ \y, - z,\^"^-^dydzdudvdsdr. (3.6) 

By Lemmas 17.41 and 17.51 we have the estimate 

p,{X^ - y)pAX^ - z) 



d 

xips{t - s - u)Lps'{t- r - v)\u- wp-^°^^TT \yi - Zi\'^^'~'^dydzdudv 

< C\s-r\'"''-^f[\Xl^X:.\'"'-\ (3.7) 

1=1 

where and in what follows C > denotes a constant independent of e and d. 



The expectation of this random variable is integrable in [0,i]^ because 
E"" f f\s-r\'"--^f{\Xl-Xl\^"'-^dsdr 



^0 



2=1 
nt nt 



where 

1 •* 
K = 2Ho + - V(2i?., - 2) - 1 > 0. (3.9) 

z=l 

and ^ is a standard symmetric a-stable random variable. 

As a consequence, taking the mathematical expectation with respect to 
X in Equation (I3.6p . letting e — e' and J = (5' and using the estimates p.7p and 
dM]) yields 



^ 2 



^t: 



<c. 

u 



This implies that almost surely ^^'^ belongs to the space "H for all £ and (5 > 0. 
Therefore, the random variables V^'^ — VF(^('^) are well defined and we have 

E^l^{ytivlf) = E^ {Ai:tAif\^ 



For any s ^ r and Xg ^ Xr, as e, e', (5 and (5' tend to zero, the left-hand side 
of the inequality ^1}, converges to \s - r\^"'>-^ nti \^'s ~ X'^\^"'~^ . There- 
fore, by dominated convergence theorem we obtain that E^ E'^ {Vfl,. V/^' ) 
converges to St as e, e', d and S' tend to zero, where 

Thus we obtain 

2 / c.\2 / c , ,,\ /,„\2 



i? (^*!f - <''' ) - i? (^t!f ) - 2i? (yt!;^^:?') + e {v^/ )' 



0. 



This implies that V/^' " is a Cauchy sequence in L^ for all sequences e„ and 
(5„ converging to zero. As a consequence, V/j' " converges in L^ to a limit 
denoted by Vt,x, which does not depend on the choice of the sequences £„ and 
Sn- Finally, by a similar argument we show p.Sp . ■ 



Proposition 3.2 Suppose 2Ho + — J2i=i{'^Hi~2) < 1. Then, conditionally to 
X the family V/^ does not converge 
a non-zero set of trajectories of X . 



X the family V/^ does not converge in probability as e and S tend to zero, for 



Proof We will prove by contradiction. Suppose ViJ^ converges to Vt^x in 
probability as e and 5 tend to zero given X for almost all trajectories of 
X. Since given X, V^i is Gaussian, so E'^iy^^) = lim £;'^(y/^'')2 = 

' ' e— >-0,i5— >0 ' 

/o /o I* - r?"°~^Ilt=i \^i - Xl\^"'-^dsdr < oo for almost all trajectories 
of X. As in the Step 2 in the proof of the Theorem 13.31 we can prove that 

(/o /o 1^ " r|2^o-2 -Qd^^ i^i _ X;|2«'-2^sdr) ' is sub-additive and it's also fi- 
nite almost surely, hence we have 

Eexp( f* f' \s - rl^^^-^HlXl - Xi.\^"^-^dsdr\ < oo, 

implying that E J^ J* \s ~ r\^"°^^ULi \^l - Xl?"'^'^dsdr < oo which is a 
contradiction. ■ 

The next result provides the exponential integrability of the random vari- 
able Vt^x defined in p.4p . 

Theorem 3.3 Suppose that 2Hq + i J2t=ii'^Hi - 2) > 1. Then, for any XeR, 
we have 

Eexp(\ I I S{X^^^ ~ y)W{dr,dy)] <oo. (3.10) 

V Jo JR-i J 

Proof The proof will be done in several steps. 
Step 1 From (|3.5p we obtain 



i?e^^*- = E^ exp (^an f^ f^ \s - r^-^ -Q |^. _ Xl\^^^-^dsd^ , 

and the scaling property of the stable Levy motion yields 

Ee^v,,^ = Ee^"^, (3.11) 

where /i — ^anf^^^, where n has been defined in p.9p . and 

Y= / \s-r\^"°-^T\\Xl-Xi.\^"'-^dsdr. (3.12) 

Jo Jo ,^1 

Then, it sufhces to show that the random variable Y has exponential moments 
of all orders. 
Step 2 Let 

Zt = (f f \s - rp«"-2 [] \Xl - Xl\^"^-'dsd)\ . 



Wc use the the identity 



\s-r\ 



2Ho-2 



2HQ-3 2HQ-3 



Co \s - u\ 2 |r - u| 2 du 



\K - X, 



i\2Hi-2 



^a \Xl-x\^^\Xl~x\^^dx i = l,---,d 



1/2 



where Ci only depends on Hi for i = 0, 1, • • • , d. 
We have 



where 



Ct(w,Xl,--- ,Xd) = 



^^(m,xi,--- ,Xd)dudxi---dxd 






a;,- 1 2 ds 



i=0 



For ^1,^2 > 0, by triangular inequality 



Zti+t2 < ^tr 



^ti+t2(w,xi,-- • ,Xd)-6i(",a;i,- •• ,Xd) dudxi ■ ■ ■ dx. 



1/2 



Write XI ^Xl^.^X^. 



iti+t2{u,Xi,--- ,Xd) -^ti{u,Xi,--- ,Xd) 



i=0 
d 



tl+t2 



l[\Xl^x,\^ds 



i=0 



t2 



s + h-u\^l[\Xl+Xl 



Xi I 2 ds 



By translation invariance, 



Cti+t2{u,xi,--- ,Xd) -6i(w,a;i,--- ,Xd) 
^? (u, xi , • • • , Xd)dudxi ■ ■ ■ dxd 



where 



i^:,{'^,Xi,--- ,Xd) = [Wc 



i=0 



2Hq-3 
Is — U\ 2 



n~. 2H:-3 

\Xl-x,\-^ds 



Therefore, the process Zt is sub-additive, which means that for any ^1,^2 > 0, 
■^ti+t2 ^ ^ti + -^t2) where Zt^ is independent of {Zs^ s < ti} and has a 
distribution same as Zt^ ■ 

Step 3 Notice Zt > is non-decreasing and path-wise continuous. By Theorem 
1.3.5 in g], for any 6* > and i > 

E exp { 9Zt } < 00 

and the hmit 

hm - log E exp \ 6 Zt] = -^ie) 

where < *(6') < 00. By the scahng Zt = f/^Zi with 

1 '^ 
j^2Ho + -y{2H,-2), (3.13) 

i—l 

we have 

hm -logEexp{0Zt] ^ lini ^log^explfi^Z 2)0^ == ^(l)^^ 

Using Chebyshev inequahty, we have e^*P|Zf > t} < Ee^^K Hence we 
have 9t + logPJZt > t) < logEe"^*, and then hmsup^^^ ^ logPJZt > i} < 

hmt_^oo jlogEe'^^* -6 = '<i'{9)-0 = 0~'^{l)-9 which is strictly negative when 
we choose 6* > sufhciently smaU. Hence there exists C > such that 

hmsup - logPJZt > t} < -C. 



i— ^00 

So we have the bound 



4Zi >t^\ <exp{-CO 



Since for random variable X>0, Ee^ ^ E J^ e^dy + 1 = J^ ¥{X > yje^dy + 
1, we have 



Ee'^^ 



pea _2_ 

< / V{0Z^-'' >y}e^dy + l 

°° ^ 

< Y^ F{0ZI-'- > K}e^+^ + 1 

K=0 

00 



K=0 



10 



This give the critical integrabihty 

Eexphz^X 



< oo 



for some 6 > 0, which implies that Eexpi \Zf> < co for all A > since 
1 <7<2. 



4 Feynman-Kac formula 

We recall that H^ is a fractional Brownian sheet on K_|_ x W^ with Hurst param- 
eters {Hq, Hi, ... , iJrf) where Hi e (^, 1) for z = 0, . . . , d. For any e, S > we 
define 

W''^{t,x):= [ [ ips{t-s)p,{x-y)W{ds,dy). (4.1) 

Jo Jw 

In order to give a notion of solution for the heat equation with fractional noise 

(jl.ip we need the following definition of the Stratonovitch integral, which is 

equivalent to that of Russo- Vallois in [T7] . 

Definition 4.1 Given a random field v — {v{t, a;), i > 0, a; G M"*} such that 

I \v{t, x)\dxdt < oo 
Jr-^ 

almost surely for all T > 0, the Stratonovitch integral L L^ v(t, x)W{dt, dx) is 
defined as the following limit in probability if it exists 

lim / / v{t,x)W^'^{t,x)dxdt. 
We are going to consider the following notion of solution for Equation 

Definition 4.2 A random field u = {u{t,x),t > 0,.t e R''} is a weak solution 
to Equation fl.l]) if for any C°° function if with compact support on W^ , we 
have 

u{t,x)tp{x)dx — / f{x)ip{x)dx — — / / u{s,x){—A)^ip{x)dxds 
' Jr'' 2 Jo jRd 

+ / / u{s,x)ip{x)W{ds,dx), 
Jo Jr'' 

almost surely, for all t > 0, where the last term is a Stratonovitch stochastic 
integral in the sense of Definition \4-l\ 

The following is the main result of this section. 

11 



Theorem 4.3 Suppose that 2Ho + — X]i=i(2-f^j^2) > 1 and that f is a bounded 
measurable function. Then process 



u{t,x) = E^ UiXnexp 
is a weak solution to Equation il-l]) . 



S{X^_,-y)W{dr,dy) 



(4.2) 



Proof Consider the approximation of the Equation (jl.ip given by the foUowing 
heat equation with a random potential 



du 



s,S 



dt 

u'''iO,x)^fix) 



{-A)'^u^-^ + u^'^Wt'^ 



(4.3) 



u'-'\t,x)=E'' (/(XHexpl / W''\t-s,X^^)ds 

where X'^ is a d-dimensional symmetric a-stable Levy motion independent of 
W starting at x. By Fubini's theorem we can write 



^From the classical Feynman-Kac formula we know that 

/ / /■* 



W'^\t-s,X^^)ds 



Mt-s- r)p,{Xf - y)W{dr, dy) ds 



V^O 
t n / pt 



Mt-s- r)p,{X^ - y)ds W{dr, dy) 



v: 



where Vf" is defined in (13.31). Therefore 



u-\t,x)^E^[f{x-)^^v{yii)) 



step 1 We will prove that for any a; € R'' and any t > 0, we have 

lim ^^|u"^*(t,a;)-w(t,a;)|P = 0, 



(4.4) 

for all p > 2, where u{t, x) is defined in (|4.2p . Notice that 
E'^\u^^\t,x) - u{t,x)\'^ = E"^ \e^ {^f{Bf) [exp {yl'^) - exp (l^,,,)]) [ 

< 11/11^ i? |exp(K!;')- exp (Ft,,) [, 

where Vt^x is defined in p. 41) . Since exp ( Vf'^i ) converges to exp [Vt^x) in prob- 
ability by Theorem 13. 11 to show (|4.4p it suffices to prove that for any A G K 

(4.5) 



sup -E exp ( Ay^ a. ) < oo. 

e,5 



12 



The estimate (g^]) follows from dSS]), dST]), and (pm)) : 



i?exp(AF,^i^)=£;exp(^- 

<Ee.p(^Cm\r-sn-'f[\X^^-X: 



2Hi-2 



drds 



< CXO. 



(4.6) 



Step 2 Now we prove that u{t, x) is a weak solution to Equation (|l.ip in the 
sense of Definition 14.21 Suppose 1^9 is a smooth function with compact support. 
We know that, 



u'^' {t,x)ip{x)dx 



1 



f{x)ip{x)dx — — / / u'^' {t,x){—A) ^ ip{x)dxds 

2 Jo jR'i 



+ 



'0 JS.'' 

Therefore, it suffices to prove that 



t r ft 

irs.S 



u^'\t,x)ip{x)W^'\s,x)dsdx. 



(4.7) 



lim / / u'^'^{s,x)(p{x)W'^'''{s,x)dsdx — I I u{s,x)(p{x)W{ds,dx), 

in probability. From (|4Jl) and (|44)) it follows that /^ J^^ u'-^is, x)ip{x)W'''^{s, x)dsdx 
converges in L^ to the random variable 



G— I u{t,x)ip{x)dx — / f{x)(f{x)dx + 



as e and S tend to zero. Hence, if 



.<^,S/ 



1 



u{t, x){—A) 2 ip[x)dxds 



ie,5. 



Be,s — I / ("u^' {s,x) — u{s,x))(p{x)W^' {s,x)dsdx 

'0 JR'i 



converges in L^ to zero, then 



m(s, x)ip{x)W^' dsdx 



u^'\s, x)ip{x)W'''dsdx - B^^s 



converges to G in L^. Thus u{s,x)Lp{x) will be Stratonovich integrable and 

u{s,x)ip{x)W{ds,dx) — G, 

JR'' 

which will complete the proof. In order to show the convergence to zero of 
B^^S: we will express the product (u^^*(s,x) — u{s,x))W^'^{s,x) as the sum of a 
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divergence integral plus a trace term (see (|2.4I) ) 



{u^'' (s, x) — u(s, x))ips [s — r)pg {x — z)5Wr, 



+(D(u^'*(s, x) - u(s, x)), ips{s - ■)pe{x - ■))-H ■ 
Then we have 



e,5 = I / (Pr:z0^^r,z 
JR'' 

+ Lp(x){D{u'''^{s,x) -u{s,x)),Lps{s- ■)p^{x- ■))udsdx 

Jo Jm'' 

= Bis + Bis, (4.8) 

where 

't'P,z — / / l""^' {s,x) — u{s,x))^{x)(ps{s ~ r)p^{x — z)dsdx, 
Jq Jr-' 

and 6{(j)'^-'^) = /^ /j^^ 4>rfz^^r,z denotes the divergence or the Skorohod integral 
of^'^-''. 

Step 3 For the term B^ g we use the following L^ estimate for the Skorohod 
integral 

EiiBlsf] < E{U'''fn) + EiWDf^'Wu^u) ■ (4-9) 

The first term in (14.91) is estimated as follows 



EiW^'Wl) - fl f f E[{u^^'{s,x)-u{s,x)){u^^'{r,y)-u{r,y))] 
Jo Jr'' Jo Jr-^ 
XLp{x)(p{y){ips{s - ■)pe{x- ■),ips{r~ ■)pe{y - ■))-Hdsdxdrdy. (4.10) 

Using Lemmas 17.41 and 17.51 we can write 

{ips{s - ■)pe{x - ■),(ps{r - ■)Pe{y - ■))n 

= an \ ips{s~ cr)(ps{r -T)\a ~ Tl'^^^^'^dadr] 

\J[a^t]^ J 



X / Pe{x~ z)ps{y~w)Y[\zt- 

\Jr2^ r=i 

< C\s~rn-'f[\x-yr^-\ (4.11) 



dzdw 



for some constant C > 0. As a consequence, the integrand on the right-hand 
side of Equation (|4.10p converges to zero as e and S tend to zero for any s, r, 
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X, y due to (|4.4I) . From (|4.6p we get 

sup sup sup E (^u'^'^{s,x)) < 11/11^ sup sup sup Eeup (2V^'^) < oo. 

E,S x<£R''0<s<t e,(5 £ceR<'0<s<t 

(4.12) 
Hence, from (I4.1ip and (|4.12p we get that the integrand on the right-hand side 
of Equation (|4.10l) is bounded by C\s — r\'^^°~'^ Yli^i \xi — 2/iP^'~^, for some 
constant C > 0. Therefore, by dominated convergence we get that £^(||0^'''||^) 
converges to zero as e and 6 tend to zero. 

Step 4 On the other hand, we have 



?x 



D{u'-'{t, x)) = E^ fiXn eMVt%")At 



^s,5\ *6,5 



where A^'^ is defined in p.2p . Therefore, 



,e'.S' , 



E{D{u'^\t,x)),D{u'^'{t,x)))n 
= E'^E^(^f{Xl+x)f{X^+x) 

xexp(K!;'(x^) + K!;'(^'))«'f(^'),</(^')>«), (4.i3) 

where X^ and X^ are two independent d-dimensional symmetric a-stable Levy 
motions, and here E^ denotes the expectation with respect to [X^^X"^). Then 
from the previous resuhs it is easy to show that 

Yim E{D{u'^\t,x)),D{u'' ''' (t,x)))H 



E[f{Xl+x)f{X^+: 



X exp 



Uh 



E 



j,k= 



f f\s^r\''"^'-''\\\Xi^'-X^^'\''"^-^dsdr\ 
1-^0 ^0 ,^1 j 



Jo Jo i^i 



(4.14) 



This imphes that u'^'^{t,x) converges in the space H)^'"^ to u{t,x) as S I and 
e I 0. Letting e' = e and S' = S in (|4.13p and using the same argument as for 
(|4.12p . we obtain 



sup sup sup _E D(u'^' {s,x))\\ <oo. 

e,S xGR''0<s<t 



Then 



E\m'''\\n^n 



E{Diu'''{s, x) - u{s, x)),D{u'''{r, y) - u{r, y)))n 
Jw^ Jo 

xip{x)(p{y){(ps(s~ ■)pe{x- ■),(ps{r- ■)pe{y - ■))ndsdxdrdy 
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converges to zero as e and S tend to zero. Hence, by (|4.9p B^ g converges to 
zero in L^ as e and S tend to zero. 

Step 5 The second summand in the right-hand side of (|4.8p can be written as 

B'^g = if{x){D{u^'\s,x)-u{s,x)),(fsis~ ■)peix- ■))-Hdsdx 

Jo JR'' 

f ^{x)E'' (fiX^) exp (1/%^) {A',i,Ms - ■)P.ix - ■))h) dsdx 
Jw^ 

[ [ ^{x)E'' {f{Xf) exp {Vs,,) {6{X'^_. - ■),Ms " ■)Pe{x " ■))«) ^sda; 
/o Js. 

33 d4 



■^e.5 - B s 



where 



J[0,s]3 jR2d ^-^ 

Xips{s — v)p^{x — z)dydyzdrdrdv , 
and 

{6{X!_.~-),Ms~-)Pe{x--))n 
= au I I 1^'^°-' n \^r' ~ y^\''"''^m{r - v)p,{x - j;)d2/di.rfr. 

Lemma 17.41 and Lemma 17.51 imply that 

{Ali, ^sis - ■)Pe{x ■))u <C r r'""-' f[ \Xl\^"^-'dr, (4.15) 

and 



X - ■))u <C r r2^»-2 fr \Xl\^"^-'dr, (4.1( 
•^0 *=i 



((5(i?f_.-.),^5(s--)Ps( 



for some constant C > 0. Then, from (|4.15p and (I4.16P and from the fact that 
the random variable /„ r^^o-'^ Y[i=i \X^\'^^^^'^dr is square integrable because of 
Lemma l7.61 we can apply the dominated convergence theorem and get that B^ g 
and B'^ g converge both in L^ to 

an f I ^ix)E^ |/(Xf)exp(y,.,) I r^""-'Y\\Xl\'"^-'dr] dsdx, 



f I ^{x)E^ [/(Xf)exp(y,,.) /V^^"-^n|X 
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as e and i5 tend to zero. Therefore B^ g converges in L^ to zero as e and S tend 
to zero. This completes the proof. ■ 

Remark. The uniqueness of the solution remains to be investigated in a 
future work. The definition of the Stratonovich integral as a limit in probability 
makes the uniqueness problem nontrivial, and it is not clear how to proceed. 
As a corollary of Theorem 14.31 we obtain the following result. 



Corollary 4.4 Suppose 2Hq + — X]j=i(2^j ~ 2) > 1. Then the solution u(t, x) 
given by (RT^ has finite moments of all orders. Moreover, for any positive 
integer p, we have 

E{u{t,xY)=Elf{f{Xl+x) (4.17) 



X exp 



^E f f\s-rr°-^Y[\x^^-x';-r^-Hsd. \ 



where Xi, . . . , Xp are independent d-dimensional standard symmetric a-stable 
Levy motions. 

5 Holder continuity of the solution 

In this section, we study the Holder continuity of the solution to Equation (|l.ip . 
The main result of this section is the following theorem. 

Theorem 5.1 Suppose that 2Hq + — J2i=ii'^^-i. — 2) > 1 and let u{t,x) be the 
solution of Equation lil.l]) . Then u{t,x) has a continuous modification such 
that for any p e (O, ^) (where k has been defined in Ii3.9\) ). and any compact 
rectangle I C K+ x R*^ there exists a positive random variable Kj such that 
almost surely, for any (s, x), {t, y) £ I we have 

\u{t, y) - u{s, x)\ < Ki{\t - s\P + \y- xD. 

Proof The proof will be done in several steps. 

Step 1 Recall that Vt.x = L J^d ^{^t-r ~ y)W{dr, dy) denotes the random 
variable introduced in p.4|) and 

u{t,x) = E''{f{Xnexp{V{t,.)). 

Set V = Vs^x and V = Vt,y. Then we can write 

E^\u{s,x) - u{t,y)\P = E^\E^{e'^ - e^)|P 
< E"^ {E'^WV -V\e^''^(^^^^]Y 
<E'^[{E''e^^''^^^'^^Y'\E''{V-VfY'^] 
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Applying Minkowski's inequality, the equivalence between the L? norm and the 
L^ norm for a Gaussian random variable, and using the exponential integrability 
property p.lOp we obtain 

E^\u{s,x)-u{t,y)\P <C[E^{E^{V-VfY]^ 

<Cp[E''E'^\V~V\'Y''\ (5.1) 

In a similar way to p.Sp we can deduce the following formula for the conditional 
variance oiV — V 

\Jo Jo i^i 

+ /* /* \r vn-2 -Q |^._^ _ Xl_X"^-'drdv 
Ja Jo i^i 

-2 r f \r - .^0-2 n \XU. XU + -. - V^r^'^drd^ 
Jo Jo ^^1 / 

■.^ aHC{s,t,x,y). (5.2) 

Step 2 Fix I < j < d. Let us estimate C{s,t,x,y) when s = t, and Xi ~ yi 
for all i ^ j. Wc can write 

C{t,t,x,y)^2 f f \r-vr^\{E{\^\^"^-^)E{\e"^-^-\z + e"^-')drdv, 
Jo Jo ,^^. 

(5.3) 

where z — — y ^^id ^ is a standard a-stablc variable. By Lemma 17.71 the 

\r — v\~ 

factor i?(|^p^J^^ ~ k + ^p^J^^) can be bounded by a constant if |r — wj < 
(a;^- —j/j)", and it can be bounded by C|a;j —j/jp|r — v|~" if |r — s| > (xj — yj)". 
In this way we obtain 

C{t,t,x,y)<C f Ir-vl'^-^drdv 

J {0<r,v<t,\r~v\<{xj-yj)''} 

+C\xj-yj\^ f Ir-vl^-^-idrdv 

J {0<r,v<t,\r-v\>(xj-yj)'^} 

<C\x,-y,r- 
So, from (|5.ip we have 

E'^\uit,x)~u{t,y)\P<C\x,-y,\^^''. (5.4) 
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Step 3 Suppose now that s < t, and x — y. Set 5 — — X]i=i(2^i ^ 2). We 
have 



C(s,i,x,x) 

ft rt 

--C 



\r — v\'^' drdv 



+ 11 \r-v\^"°-^{\r-v\^ -\r-v + t-s\^)drdv 
Ja JO 

The first integral is 0{{t — s)'^+^), when i — s is smaU. For the second integral 
we use the change of variable a = v — r,v = t, and we have 

/•s rt 

\r - w|2^o-2^|^ _yf _\r-v + t- s\^)drdv 



-Jo'^JJ^ 



^"°-^\\<7\'-\a + t-s\'\da 
Jo 



=t 



+ i-af""-^\i-a)^ -i<7 + t-sY\da 

Js-t 

=t[A' + B' + C']. 
For the first term in the above decomposition we can write 



,2Ho-2(^S 



{a^ -{a + lY)da 



'o 

/•OO 

Jo 



because 2Hq — 2 + 5— 1<— 1. Similarly we can get that 



At last, 



So we have 



/OO 



C < \ a^"°-^{a^ + {t-s-a)^)da = C{t-s)^. 
Jo 

E^\u{s, x) - u{t, y)\P < C{t - s)iP. 



(5.5) 



Step 4 Combining Equation 15.41 and Equation 15.51 with the estimates (|5.ip 
and (|5.2p . the result of this theorem now can be concluded from Theorem 1.4.1 
in Kunita [13] if we choose p large enough. ■ 
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6 Skorohod type equations and chaos expansion 

In this section we consider the foUowing heat equation on M'^ 

u{0,x) = fix). 

The difference between the above equation and Equation (|l.ip is that here we 
use the Wick product o (see [H], for example). This equation is studied in (TU] 
in the case Hi = ■ ■ ■ = Hd = 1/2. As in that paper, we can define the following 
notion of mild solution. 

Definition 6.1 An adapted random field u = {u{t^ a;), i > 0, a; G M''} such that 
E(v?'{t,x)) < oo for all {t,x) is a mild solution to Equation i6.6\) if for any 
{t,x) £ [0,oo) X E'^, the process {qt-s{x — y)u{s,y)l[Qt]{s), s > 0,y G M.'^} is 
Skorohod integrable, and the following equation holds 

u{t,x) ^ qtf{x)+ / qt-s{x~y)u{s,y)SWs,y, (6.7) 

Jo JR'^ 

where qt{x) denotes the density function of Xt and qtf{x) — j^^ qt{x — y)f{y)dy. 

As in the paper [TI^ the mild solution u{t, x) to (|6.6p admits the following 
Wiener chaos expansion 

oo 

u(t, a;) = ^ /„(/„(•, t, a:)), (6.8) 

n=0 

where /„ denotes the multiple stochastic integral with respect to W and /„(•, t, x) 
is a symmetric element in "H*^", defined explicitly as 

/„(si,yi,...,s„,2/„,i,x) = — (6.9) 

n\ 

X9t-5„(„)(a; - 2/^(n)) • ■ ■ qs„f2)-s,^t) {ya{2) " y<T(l) )gs,(i) /(2/o-(l) ) ■ 

In the above equation a denotes a permutation of {1, 2, . . . , n} such that < 
■So-(i) < • • • < s^(„) < t. Moreover, the solution if it exist, it will be unique 
because the kernels in the Wiener chaos expansion are uniquely determined. 
The following is the main result of this section. 

Theorem 6.2 Suppose that 2Ho + — X]i=i(2-^j^2) > 1 and that f is a bounded 
measurable function. Then the process 



u{t,x) = E 



X 



f{Xn exp I / * / S{X^_, - y)W{dr, dy) 



t-r 




Ian f f \r - s^^^^'-^fWxi- Xl\'"^-' drds 



2 



Jo 



i=l 

is the unique mild solution to Equation lil.l\) . 
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(6.10) 



Proof /^From Theorem 13.31 we obtain that the expectation E^ in Equation 
(|6.10p is well defined. Then, it suffices to show that the random variable u{t, x) 
has the Wiener chaos expansion (j6.8p . This can be easily proved by expanding 
the exponential and then taken the expectation with respect to X. 

Theorem 13.11 implies that almost surely 5{Xf_. — •) is and element of 
T-L with a norm given by (|3.4[) . As a consequence, almost surely with respect 
to the stable Levy motion X, we have the following chaos expansion for the 
exponential factor in Equation (j6.10p 

expf / / 5{Xf_,-y)W{dr,dy) 

T rt rt d \ CO 

^lan / \r-sn-^Yl\x:.-Xl\'"'-'drds]^J2lr.ign), 

where gn is the symmetric element in "H®" given by 

gn{si,yi,...,Sn,yn,t,x) = — ^(Xf_ - yi) ■ ■ ■ d{Xf_ ~yn)- (6.11) 
Thus the right hand side of (|6.10p admits the following chaos expansion 

OO -. 

u{t,x) = Y^—In{K{-,t,x)), (6.12) 

with 

K{t, x) = E"" [f{xr)5{X-^.s, - 2/i) • • • 5{X't_,^ - y„)] . (6.13) 

This can be regarded as a Feynman-Kac formula for the coefficients of chaos 
expansion of the solution of (|6.6|) . To compute the above expectation we shall 
use the following 

E"" [f{Xf)5{X't - 2/)|-F,] - j ^ qt-siX't - z)f{z)5{z - y)dz 

= qt-s{X't-y)f{y). (6.14) 

Assume that < Scr(i) < • • • < Scr{n) < t for some permutation cr of {1, 2, • • • , n}. 
Then conditioning with respect to -^t-s„(i) and using the Markov property of 
the Levy motion we have 

K{t,x) = ii;^{i?^[<5(xf_ -y.(„)) 



x---5(^t%„,,, "2/.(i))/(^r)|-^*-.„<j} 



E 



(1 

X 



'^(^t-..(„, - y<n)) ■ --Hx^-s^,,, - y.iDhs^aJi^t-s^,,,) 
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Conditioning with respect to J^t-s^,2) ^^"^ using (J6.14I) . we have 



= E- KXt-s^^^^ - y.(„)) • • • 5{Xts^^,^ - 2/.(2)) 



xE 



X 



'5(^.-.„<„-jy.(i))9..a,/(^r-.„<„)|-^*- 



S^(2) 



= E 



X 



^i^t-sl^^^ ~ y<r(n)) • • • S{X^_,^^^^ - 2/^(2)) 



Continuing this way we shall find out that 

h„{t, x) = qt-s^f^-, {x - y<,(„)) • • • <?s„(2,-s„(i, (y<T(2) - ya{l))qs,^,■, /(y<T(i)) 

which is the same as (|6.9p . ■ 

Remark. The method of this section can be applied to obtain a Feynman-Kac 

formula for the coefficients of the chaos expansion of the solution to Equation 



with 
K{t,x) 



E 



X 



OO -J 

U{t,x) = V —In{hn{-,t,x)) , 
n=0 

f{Xt)5{X't_,^^y,)---5{Xts„-yn) 
e.^(\aHff\r-sr''-'i{\Xl-Xl 



2Hi-2 



drds 



(6.15) 



From the Feyman-Kac formula we can derive the following formula for 
the moments of the solution analogous to (|4.17p . 



E{u{t,x)n = E( i[f{X{+x) 



^j=i 



X exp 



an 



E 



j,k=lj<k 



^0 



\2Hn-2 



Y[\Xi^' - x!;-'\^"'-^dsdr 



where p > 1 is an integer, and X^ , 1 < J < rf, are independent d-dimensional 
stable Levy motions. 
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7 Appendix 

Lemma 7.1 For any deterministic sub-additive function a{t), t G M+. the equal- 
ity 

lim t^ a(t) ~ inf s^ a(s) 

t-J-oo s>0 

holds in the extended real line [—00,00). 

Lemma 7.2 Suppose < /3 < 1, e > 0, a; > 0, and that X is a standard normal 
random variable. Then there is a constant C independent of x and e (it may 
depend on j3) such that 

E\x + eX\-'^ < Cmm{e-'^,x-^). 
Proof It is straightforward to check that K ~ sup2>o -^k + X\^^ < co- Thus 
E\x + eX\-^ = e-^E\- + X\-l^ < Ke~^ . (7.16) 

On the other hand, 

Elx + eX]-/^ ^ -^ f Ix + eyl'^e-'^dy 
V27r Jr 

Air 2 

X + ey\ ^e 2 dy 



V2vr \J{\x+ey\>^} 

+ \x + ey\''^e-'^dy] . 

J{\x+ey\<^} J 

It is easy to see that the first integral is bounded by Ca;"^ for some constant 
C. The second integral, denoted by B is bounded as folfows. 

B = C- \zr^e-^^dz<C- \z\-^e^^dz 



= C-e-^x-P <Cx-^ . 
e 

Thus we have E\x + €X\-^ < C\x\~^ . Combining this with ((7?TB| . we obtain 
the lemma. ■ 

Lemma 7.3 Suppose < /3 < 1, e > 0, a > 0, anii t/iat K is a standard symmet- 
ric a-stable distributed random variable. Then there is a constant C independent 
of X and e (it may depend on /3) such that 

Elx + eYl-l^ <Ce'^. 
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Proof 



E\x + eY\-'^ 



= /-^{k-r^}(Oe"«e-l«l°dC 

■/R 

i-F{|-r^}(^)e"«e-l«l°ci^ 
i|i|/3-ie-«e-l«l°de 



Lemma 7.4 Suppose a € (0, 1). There exists a constant C > , such that 

sup / p^{xi + yi)pe'ix2 + y2)\yi ~ y2\~°'dyidy2 < C\xi - X2\~°'. 
Proof We can write 
/ Peixi+yi)pe'ix2 + y2)\yi-y2\^"dyidy2^ E (\eXi-xi-e'X2 + X2\^°') ■ 

JR2 

Thus Lemma [7^ follows directly from Lemma [7.21 ■ 

Lemma 7.5 Suppose a G (0, 1). There exists a constant C > , such that 

sup/ / ipsit- si-ri)tps'{t- S2-r2)\ri~r2\~"'dridr2 <C\si- S2\~°' 
S.S' Jo Jo 

Proof Since 

1 _xi 1 

PS{X) >psix)I^^ys^ix) = -7==e 2^/[oy5](x) > -j==ip^{x) , 



/2ttS '■^"^ ~ V2TTe 

the lemma follows from Lemma 17.41 



Lemma 7.6 Suppose that 2Hq + i J^t^iC^^i - 2) > 1. Let X^,...,X'^ be 
independent one- dimensional symmetric a-stable Levy motion. Then we have 



Elf\'""~'f[\Xl\'"'-'ds\ <oo. 
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t rs 

■2Ha-2 



Proof Wc can write 

Cd \ -^ 



d 
X 

1=1 



Let y be a standard symmetric a-stable distributed random variable. ^From 
Lemma l7.3l taking into account that 2 — 2Hi < 1, we have when r < s, 

< CE[\Xl\^"^-^{s-rf-^^) 

2H.:-1 2H,-2 

< Cr^^{s-r)^^. (7.17) 
As a consequence, the conclusion of the lemma follows from the fact that 

Jo Jq 
because 2i7o - 2 + i Eti(2^« - 2) > -1 and i Eti(2^« - 2) > -1. ■ 
Lemma 7.7 For any Q < /3 < 1, 

for some constant C > 0, where y > and ^ is a standard symmetric a-stable 
random variable. 

Proof Notice first that £'(|^|"'^ - \y + ^\^^) < C where C > is a constant, 
since limy_j.oo E\y + CI " = 0. 
On the other hand, 

E{\^\-^-\y + i\-P) 

=C f \Cf-^e-\i\"il-e'y^)d^ 

Jr 

=c/|Crie-l«l°(l-cos(yOK 

JR 

<C f \if-'e-\i\"y'ed^ 

Jr. 

<Cy^ 
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